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Abstract 


Ws Investigate tha feasibility of producing the observed velocities 
and teaperatures of solar wind heavy ions (v^ * ^p ^A* * ^^p^ 

resonant cyclotron interaction with Icf t-polarised hydroaagnetic waves. 

We set up a favorable case" scenario. The waves are parallel-pro- 

pagating and dispersionless. The energy for the wave acceleration and 
heating is taken from saturated low-frequency Alfvra waves via a cascade 
to higher frequencies. This scenario is incorporated into a numerical 
solar wind code and agreement with observation is tested. We find that 
the resonant cyclotron interaction fails on at least three points, even 
in this most favorable case. First, we show that no preferential heating 
or acceleration is obtained for heavy ions unless the wave power spectrum 
is steeper than which is not often observed. Secondly, we find that, 
when sufficiently steep spectra are assumed, the resonant cyclotron inter- 
action is unable to produce helium or oxygen velocities in excess of a 
fraction of the Alfven speed with respect to the protons inside 1 AU. The 
simultaneous wave Interaction with the protons prevents the heavy ions from 
attaining flow speeds of v^ 4- and this effect is understood in the form of 

a limiting differential apeed. This limiting speed depends on A/Z of the 
heavy ion. Solar wind iron ions, with larger A/Z than helium or oxygen, 

can reach v 4- V. in this model. Thirdly, the preferential wave heating is 
P A 

insufficient to produce T^ • ^j,Tp at 1 AU, and we find th»»t T^^/T^ is nega- 
tively correlated with proton speed, the opposite of the observed behavior. 




I. Introduction 

The properties of solar wind ions heavier than H have presented a 
pussle to theorists for some time (see Heugebauer » 1981 Stb; HolIwet « 

1981a for recent reviews). The observation that Che heavy ion velocity » 
frequently exceeds the proton velocity, v^, in fast solar wind 
( Formisano. et a l. . 1970; Bollea, et al .. 1971; Asbridge, at al .. 197A, 

1976; Hirshberg, et al ., 1974; Ogllvle , 1975; Bosqjed, et al .. 1976; 

Marsch. et al .. 1982a; Ogllvie, et al ., 1982) has not been satisfactorily 
explained. Equally pussling is the fact that the heavy ion temperatures 
are also larger than that of the protons, with the temperatures often 
appearing proportional to the masses of the ions, l.e. ^ where 

is the atomic weight of the heavy ion (Robbins . et al . , 1970; Ogllvie 
and Zwally . 1972; Bosgued, et al .. 1976; Ogllvie et al .. 1980, 1982; 

Schmidt et al . , 1980; Bochsler and Geiss . 1982) 

The most abundant heavy ion. He , has been studied most extensively 
and positive correlations have been found between the differential speed 
of the alpha particles with respect to the protons, v^ - v^, and both the 
proton and Alfven speeds (Asbridge et al .. 1976; Neugebauer . 1981a; Grunwaldt 
and Rosenbauer . 1978; Marsch et al . . 1982a). The differential speed is also 
seen to be positively correlated with the alpha-proton temperature ratio 
(Neugebauer and Feldman . 1979). 

The correlation with the Alfven speed, V^, is usually interpreted to 
Indicate that the mechanism producing the preferential acceleration is a 
wave-particle interaction with electromagnetic waves having phase speeds 
of V^. This interpretation is strengthened by observations closer to the 


sun. Marsch et al . (1982a) report that inside 0.8AU, Che alpha particles 


are often seen to flov at the Alfven velocity with respect to the 
protons, proaptlng the laage of the helium Ions "surfing” on the waves. 
Furthermore, Isenbere and Hollwea (1982) have shown that preferential 
acceleration produced by wave dissipation Is accompanied by preferential 
heating, Implying a natural explanation of the relation bet%reen v^ - v^ 

and T /T . 

o p 

Host recent theoretical efforts to explain the observations have 
concentrated on the resonant cyclotron Interaction with parallel-propa- 
gating Ion-cyclotron or high-frequency Alfven waves. The low-frequency 
analogs of these waves are commonly observed In the solar wind (Belcher 
and Davis , 1971; Behannon . 1978). The Interaction Is known to be 
capable of accelerating heavy Ions from Vj^ < v^ to Vj^ > v^ (Hollweg and 
Turner . 1978; Dusenbery and Hollweg . 1981, hereafter deslgnatnl DH; 

Hollweg . 1981a). Additionally, the fact that the proton temperature 
In fast solar wind Is enhanced In the perpendicular direction (Bame, et al . , 
1975; Marsch. et al . . 1982b) Is conjibtent with the expected result of 
resonant ion-cyclotron heating. 

However, the efficacy of this mechanism as the explanation of the 
observed preferential heating and acceleration has yet to be determined. 
Hollweg and Turner (1978) Investigated the resonant acceleration of a 
test population of alpha particles by nondlsperstve Alfven waves. The 
helium was assumed to move In a background electron-proton solar wind 
which was unaffected by the resonant interaction or by the presence of 
the helium. Furthermore, the temperatures of all the species were pre- 
determined, so wave heating was not considered. Under these assumptions, 
the authors were able to obtain v “ v + V, for a variety of models. 
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DK reported on a parametric study which calculated both heating and 
acceleration rates for the resonant interaction with both dispersive 
and nondispersive ion-cyclotron waves. The heavy Ions were again treated 
as test particles which did not affect the wave field, and there was 
no attempt to Incorporate the resonant effects into a solar wiiui model. 

This work indicated strong acceleration and heating of heavy ions by 
the resonant interaction, but found that there could be limits to the 
local acceleration. For waves satisfying a cold plasma ion-cyclotron 
dispersion relation, the sign of the ion acceleration became negative 
when v^ - v^ reached a certain fraction of depending on the temper- 
ature (see also McKenzie and Marsch . 1982). A similar limit «ms found 
for nondispersive waves interacting with warm anisotropic ion distri- 
butions where (the subscripts &nd ' n * indicate directions 

perpendicular to and parallel to the ambient magnetic field direction, 
respectively). 

Another difficulty with the resonant mechanism pointed out by DH 

concerns the energy requirement. When alpha particles are observed to 

flow at v_ ■ V V. in the inner solar system (where V. is as high as 
Cl p A A 

150 km/sec), the kinetic energy flux of the population is a substantial 
fraction of the total flux of solar wind energy. To accelerate the 
observed numbers of alpha particles to the observed velocities requires 
the energy equivalent to at least 20Z of the observed wave field. Further, 
the resonant ion cyclotron waves only constitute the high frequency 
portion of the wave spectrum while the bulk of the power is found at 
much lower frequencies. Thus, the observed wave power in the resonant 
regime is generally several orders of magnitude below what is needed. 
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So, the Importance of this mechanism remains an open question. 

The energy problem can be circumvented, at least in principle, by in- 
voking a nonlinear cascade of energy from low frequencies (as suggested 
by DH), but other questions must be dealt with. One is the simultaneous 
resonant interaction with the protons. It is necessary to show that 
the resonant waves produce preferential acceleration and heating of 
heavy ions ;«lth respect to protons. This is a point that neither of 
the above papers considered. Then, the local resonant interaction must 
be Incorporated into £ solar wind model to gauge the effects of the 
global inhomogeneities of the system. For instance, the Alfven velocity 
in the solar wind decreases with increasing radius and a mechanism which 
produced ^ " solar radius) might be expected 

to exhibit v v + V, at .3 AU if no other tcrces are present. There 
a p A 

are, however, counteracting forces such as Coulomb collisions and non- 
resonant wave interactions, which tend to equalize the velocities of 
different populations. The relative effects of the various forces mcst 
be determined from a model calculation. 

In this paper, we investigate a scenario which deals with these 
questions. This sosnario should be considered as a "most favorable case" 
picture of the resonant acceleration of heavy ions by ion cyclotron waves 
In this picture, the energy for the resonant acceleration is siphoned 
from the low-frequency part of the wave spectrum. We will hypothesize 
a saturation of the Aifvenic oscillations, i.e. a nonlinear dissipation 
which limits the power contained in the low frequency part of the wave 
spectrum. The saturation process is assumed to feed energy to higher 


s 


frequencies by an unspecified non-linear cascade providing sufficient 
energy in the resonant range to accelerate the heavy Ions. 

In view of the fmidamental limits inherent in resonant acceleration 
by dispersive ion-cyclotron waves, we consider the interaction with dis- 
persionless, parallel-propagating left-polarized waves - essentially 
high-frequency Alfven waves. Where DH showed that cold dispersive waves 
cannot resonate with helium moving faster than <■0.2 with respect to 
the protdns, the Alfven waves can affect all particles with v^- v^ < V^, 

so the most favorable resonant process is an Interaction with high- 
* 

frequency Alfven waves. 

Finally, these wave processes are incorporated into the three-fluid, 
supersonic, spherlcally-symnetrlc, time-independent solar wind equations 
which include the radial dependence of the wave phase speed as well as 
the effects of non-resonant wave forces. Coulomb collisions, solar rota- 
tion, and gravity. The heavy ions are not treated as test particles in 
this calculation. This is particularly Important for helium since Hellos 
measurements (Marsch, et al .. 1982a) have shown that the helium population 
can carry a substantial fraction of the energy of the solar wind in the 
inner solar system. Consequently, we allow the heavy ion fluid to carry 
momeni.ur< and energy and to exchange these quantities with the other 
fluids and with the waves. 

The numerical calculations show that the resonant cyclotron inter- 
action can indeed produce some preferential acceleration and heating of 
the helium and heavier ionmand some tendency for mass-dependent heating. 
However, even for this "most favorable” case, the resonant ion-cyclotron 

interaction is Incapable of accelerating the helium population to v >v -^V, 

a p A 

at ,3AU. Furthermore, for all but the most extreme values of the parameters, 


this mechanism cannot produce v 2 v V./2 for the entire range 

<1 P A 

between lOR and lAU. Similar results are obtained for oxygen ions 
8 

though the mechanism is not as constrained by the observations in 
these cases. Numerical results for iron ions are also discussed. 

The next section of this paper contains a detailed discussion 
of the resonant wave-particle interaction and the derivation of the 
relevant equations for our "most favorable case". In Section 3» the 
computational model is described. Section 4 consists of a discussion 
of the computational results and their interpretation. Section 5 
summarizes our major points and conclusions. 


II Wave-Particle Interaction for High-Frequency Alfvep Waves 
a) “Most favorable case" scenario. 

With the Intention of Investigating the most favorable situation 
for resonant preferential acceleration and heating of heavy Ions by 
left-polarlxed waves, we consider the interaction with nondispersive 
parallel - propagating waves. Essentially, these are high-frequency 
Alfven waves propagating along the average magnetic field with a phase 
speed which is independent of the wave parameters. 

The condition for cyclotron resonance between a particle and a 
circularly polarized parallel-propagating wave is that the Doppler- 
shifted wave frequency In the frame of the particle's guiding center 
is equal to the gyrofrequency of the particle. This is 

“ - hi 


where o) and k are the angular frequency and the wavevector, respectively 
of the wave, gyrofrequency of the particle and 

the superscript 'r' denotes the resonant velocity. In the solar wind, 
we will consider both waves and particles to be moving away from the 
sun, so and v^ || are of the same sign. When the particles arc posi* 
tively charged, the upper sign in (1) is appropriate for left-polarized 
Alfven or lon-cyclotron waves and the lower sign corresponds to the 
right-polarized magnetosonic or whistler mode. Equation (1) can be 
solved for the parallel velocity of a resonant particle In terms of the 


phase speed (o/k„ : 


+ 0 ^/«) 


( 2 ) 


It is clear that a particle can only resonate with the right-hand mode 
if It is already moving faster than the phase speed of the wave. Since 
we are concerned with accelerating particles from speeds below the Alfven 
speed up to this speed, we will concentrate on the left-polarized ion- 
cyclotron waves from this point, taking the upper sign in (1) and (2). 
Possible consequences of the right-hand interaction will be discussed 
in Section 4. 

Dispersive ion-cyclotron waves have a phase speed which goes to 
zero in the plasma frame ar> m in a plasma dominated by protons. 

Ic can be seer; that the resonant velocity (2) for a heavy ion then has 

z 

a maximum which occurs when the resonant v^ve frequency is in the range 
< u < Dp. This maximum velocity is always less than V^. For an 
alpha particle in a cold proton plasma, this maximum is v ^ >0.2 V. 

for parallel waves (DH) and is only slightly larger for perpendicular 
propagation (McKenzie and Marsch , 1982). This means that a particle 
moving along the field at a speed faster than this maximum cannot resonate 
with dispersive Ion-cyclotron waves and the resonant Interaction with 
these waves cannot accelerate heavy ions to V^. 

Nondlspcrsive waves can interact with faster particles. Since we 
do not know the correct dispersion relation (which depends on the dis- 
tribution functions of all the plasma constituents), we will consider 
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the nose favorable situation and asauna tha interacting wavaa to be 
parallel-propagating Alfven %iave8, with a phase speed independent of 
wave paraoeters. There is some theoretical evidence (Gary, at al .i 1975) 
that ion-cyclotron waves remain essentially dispersionless in a warm 
plasma, but our primary motivation in using nondlspersive waves is that 
they seem capable, in principle, of accelerating heavy ions up to 

’h • + ''a- 

Another aspect of this mechanism to consider is the energy available 

in the resonant waves for the acceleration and heating of the alpha 

particles. As stated by DH, one needs about 201 of the total energy 

in the observed Alfvenic fluctuations to accelerate a number density of 

helium, n^ " C.u4n^, to v^ » v + V. at .3AU in the solar wind. Very 
o p a p A 

little of the total power in the fluctuations is present at resonant 
frequencies and there is a serious problem in obtaining any substantial 
effect on heavy ions in the solar wind if one is restricted to using 
only the wave power observed in this range (see also Schwarts, et al . , 1981) . 

There exists sufficient power in the low-frequency portion of the 
wave spectrum, however, and we propose a scenario which allows the resonant 
Interaction to utilize this energy. We point out first that the low 
frequency waves are observed to dissipate as they propagate through the 
solar wind (Villante . 1980; Bavassano, et al .. 1982a, b) though the 
mechanism is not understood. For the purpose of this investigation, we 
assume that a non-linear dissipation process saturates the wave activity 
at some relative amplitude, ■ constant, where <5|^> is the 

total magnetic variance (integrated over the entire spectrum), and B 

o 

is the magnitude of the average magnetic field. The energy lost to the 


wave field eseencially comes from Che large amplicude, low frequency 
end of the wave spectrum* but we further assume the operation of a 
cascade process which trarispcrts the energy from the saturated low- 
frequency range of the spectrum to the resonant range where it is 
picked up by the plasma through the resonant cyclotron interaction. 

The details of the resonant interaction will be given by Che quasi- 
linear expressions. The most favorable situation for the resonant inter- 
action hypothesis would make the maximum amount of energy available to 
the resonant particles. Consequently* we assume that the resonant in- 
teraction is the only means of transferring the dissipated energy to 
the plasma and that the cascade is efficient enough to funnel all the 
power lost from the low-frequency waves into the resonant fluctuations. 

The saturation and cascade processes will not be specified. We 
merely note that recent observations of fluctuation amplitudes suggest 
that the waves are saturated (Villante* 1980; Bavassano* et al .. 1982a*b). 
Moreover, the typical power law character of observed wave power spectra 
has often suggested the operation of a non-linear "turbulent" cascade in 
wavenumber (Coleman . 1968; Barnes . 1981; Hatthaeus and Goldstein . 1982). 

So we feel that the scenario outlined above is a plausible one. Finally, 
this saturation scheme (without the necessity of imposing a cascade hypo- 
thesis) has been used successfully by Hollweg (1978a) to model the high- 
speed flow of the solar wind. 

b) Macroscopic equations 

This scenario allows us to utilise the macroscopic results of Isenbcrg 
and Hollweg (1982) for multi-species acceleration and heating by Alfven 
waves. In a time-independent system where the bulk fluid velocities of 
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all species are directed along the average magnetic field (equivalent 
to the corotating fraoie in the solar wind), assumlu^ all particle species 
to be thermally isotropic, they showed that the total plasma heating 
due to wave dissipation is related to the divergence of the wave action 
flux 


V • (diS) - - E (w/kjl) ^ 


(3) 


The flux of wave action for nondispersive Altven waves was found to be 


S 


<6B2> 

4it 


k Zy f 



(4) 


^ th 

where b = . The volumetric heating rate of the 1 species is 

Q. and Its bulk velocity is v The Alfven speed is V,= B (Airlp.) 

i ^01 ^ A O £ 1 

and the center of mass speed is s where p is the density 

and the summations are over all species. The phase speed of these waves 
in a differentially flowing multi-ion plasma is 


ii)/k 


II CM 


|_A CM 1 1 oi i y 


(5) 


These results are valid for WKB Alfven waves of arbitrary amplitude and 
arbitrary angle between k and In the solar wind, most waves propagate 
out from the sun. In the same direction as the plasma, so we will take 
the upper sign in (5). 
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Note that the wave heating in a moving plasma is not given simply 
by the rate at which energy is lost by the waves. In general, the 
dissipation of waves adds momentum to the plasma as well as energy, so 
some of the wave energy must go into doing work. The conservation of 
wave action includes the effect of this additional energy loss. 

The work done by these Alfven waves in a steady solar wind can be 
related to the heating rate by the First Law of Thermodynamics. It is 
found, then, that the acceleration of a given plasma species a^, due to 
the waves satisfies 


'oi 


h 


>oi 


<5B2 


2B 2 
o 


(k 


2 



Si ^0^11 

Pi tuA,, - 


( 6 ) 


In a spherically symmetric system, where all gradients are in the radial 
direction, (6) gives the change in flow energy of an ion species parallel 
to the field as - 



+ Qi/»i 




(7) 


phr 


where the subscript 'r' denotes a radial component and V 


(m/k) (B/S ). 
II r o 
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The first tent on the right-hand-side of (7) is the non-dissipative 
wave pressure (Hollweg . 1974a, b, 1978b; Chang and Hollweg . 1976; 

Goodrich , 1978; McKensie, et al .. 1979) generalised to finite amplitude 
Alfven waves with arbitrary directio 3 of k and This pressure force 

is an important source of solar wind acceleration. However, It acts 
on differentially flowing particle populations so as to equalise 
their bulk velocities, and therefore works against preferential accel- 
eration. The second term in (7) is the acceleration which necessarily 
accompanies dissipative wave heating, and it is this term which will 
produce preferential acceleration resulting from the resonant Interaction. 

The scenario for the wave-particle interaction then proceeds as 
follows: for wave amplitudes below the saturation level, there is 

no dissipation, the heating rates are sero, and the waves evolve accord- 
ing to the conservation of wave action in a time-independent system 


7 


(w 5) 


0 


(8 


In the expanding solar wind, (8) requires that increase with 

^ o 

increasing radius inside about lAU. At some point, <6B^>/B^ will reach 

' o 

the designated saturation value and initiate the non-linear dissipation 
envisioned above. Beyond this radius, the wave amplitudes are given 
by the saturation condition and (3) determines the total wave heating. 

To proceed further, we must specify how this total wave heating is 
divided among the ion species, and this is determined by the microphysics 
of the resonant interaction. 
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c) Microscopic equations 

The quasl-llnear theory of the wave-particle interaction provides 
an expression for the rate of change of a particle distribution function 
f^(y) due to the resonant interaction with a wave power spectrum P(k) 
(e.g. DH). The heating rate is then related to the second moment of 
3f^/3t: 


Qj . 1/2 p J dy [(V|| - 8fj/3t 


For a bi-Maxwelllan distribution interacting with a power-law 
spectrum of parallel-propagating left-polarized waves, the heatings 
rate is 





“i 

r > 

- V ^ 

2 

2k T 

Hi 

k,| oi„ k„ 

^ J 



where the power spectrum of magnetic fluctuations is taken to be 
P(k) “ P^(k^/k^) the temperature anisotropy is T ^ is 

Boltzmann ^s constant. 

Dispersionless waves have a phase speed (5) independent of 
wave iraber so some terms can be factored out of the integral. The 
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heating rate can be rewritten 

Q. 


f 1 


H ] 

|2<T , 
1 

Pi 

m^c 

k 4 


y . y+2 -y 

I P k (w/k„ - V .j n 
m^c I o o ' H oir i 

. J* dx |7i -r ) X + ^ (11) 


where x 


Oi 


0) 

k " ^oi. 


1 -1 Oj’ 

‘=1 ■ 


fa) 

k. 


- V 


oi. 


If we assume that the Integral Is dominated by the resonance at 
X ■ 1 (i*e. C^»l), we find the proportionality relation 

Qi « Pi (Ai/Zi>^-2 (o)/k|| - \l + (Y-1) (Y-2D/(4Ci) + ...] (12) 

where is the ionic charge. In the numerical investigation to follow, 
we will assume that is small enough so that it is sufficient to 
keep only the leading term in (12). This is valid for the input condi- 
tions we will impose and it is likely to remain valid since the resonant 
interaction with ion cyclotron waves preferentially raises the tem- 
peratures perpendicular to the magnetic field. The calculations in 
DH showed that the parallel temperatures are even lowered in many cases 
by this interaction. Furthermore, their calculations indicate that when 

V . - V < 0.5 V. (and ■ 1.5), the heating rate is essentially in- 
ch op * A 

dependent of temperature up to 2<T||j^/mj^ “ Since the leading 

term in (12) does not depend on the temperature anisotropy, the effects 

of an anisotropic particle distribution will be neglected in the analysis 


to follow. This approximation and its implications will be discussed 
further in Section 4. 

One final point remains to be made about this scenario for prefer- 
ential acceleration and heating of heavy ions. If these preferential 
effects are to be produced solely by the resonant interaction with left- 
polarized high-frequency Alfven waves, it is necessary that * ^p^*^p* 

where is given by the leading term in (12). For is 

clear that the acceleration and heating process will only act prefer- 
entially on the heavy ions, > 1* when the wave spectrum is steeper 

than k**^. Unfortunately, observed wave spectra do not appear to be steep 
enough (Behannon , 1978; Bavassano, et al ., 1982a). Typical spectral 

indices for u<n (in the plasma frame) are on the order of y^'l.S, with 
P 

occasional observations as high as y“2 but no higher. The spectra do 
steepen to at high frequencies, but it is the range (in the 

plasma frame) which determines the heating for the case of similar ion 
speeds. 

In view of this discrepancy, we Include one more assumption in the 
"most favorable case" picture: that y>2. This is admittedly a questionable 
assumption. There are no observations at 20 Rg where the bulk of the 
interaction will take place in our model, but indications from radio 
sclntellatlon measurements are that the spectra flatten as one approachs 
the sun, rather than steepening as required ( Woo , 1981). Still, the re- 
sonant ion-cyclotron mechanism for preferential heating will not work 
without this assumption, so the "most favorable case" must incorporate it. 

The scenario to be considered is now complete. When the waves have 
reached saturation amplitude, the total plasma heating is given by (3). 
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The relative wave heating of each species is determined by the leading 
term in (12) so the heating rate of each species can be found. Finally, 
the bulk acceleration of each species due to the waves is specified by 
(7). This wave interaction will be Incorporated into the numerical 


solar wind model detailed in the next section. 


Ill Computational Model 

a) Assumptions 

To test the effectiveness of our "most favorable case" heavy Ion 
acceleration by high-frequency Alfven waves, we Incorporate the interaction 
scenario into a numerical solar wind calculation which uses the following 
assumptions: 

1) The wave-averaged system is time-independent and spherically 
symmetric. Therefore, all bulk quantities have 3/3t« 0 and vary only in 
the radial direction. 

2) The calculations will be performed in the heliospheric equator and 
in the frame rotating with the sun. In this frame, all fluid velocities 
are parallw. ^ ..ne average magnetic field and the only large-scale electric 
field is the usual radial charge-separation field. 

3) Since we are not concerned here with the origin of the solar wind 
per se , we will start the calculation at 10 solar radii (Rg) from the center 
of the sun, which is outside the sonic critical points for all parameters 

of interest. We will require that the parameter values at 10 Rg be such 
as to give reasonable agreement with observations given the simplifications 
made. This requirement allows us to choose proton input values which are 
consistent with those suggested by the wave-driven high-speed stream model 
of Hollweg (1978a). 

4) The calculation will deal with three fluid species: electrons, 
protons, and a heavy ion species. We assume isotropic distributions for all 
species, characterizing the bulk properties by a flow velocity, y^, and a 
temperature T. In the explicit calculation of the quasl-linear wave heating 
rate and of the Coulomb collision terms, the distributions are taken to be 
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Maxwelllans, but this stricter assumption Is not needed elsewhere. 

Considerable particle anisotropies are observed In the solar wind, but 
the assumption of isotropy Is in keeping with our "most favorable case" 
intent, at least as concerns the ions. Ion anisotropies will most seriously 
affect the wave-particle interaction, but the Ion-cyclotron interaction will 
tend to produce T|| which reduces the heating (and acceleration), as can 
be seen from (12). Reduced heating as a consequence of Tj_> T| Is also seen 
in DH. They found that the interaction of a pancake distribution with left- 
polarized waves results in deceleration and cooling for much of the para- 
meter space. The effect of particle anisotropies on the resonant acceleration 
and heating will be discussed further in the next section. 

5) The heavy Ion species is allowed to affect the flow; that Is, they 
are not test particles. They interact with the other species by Coulomb 
collisions, participate in determining the spiral field angle, the phase 
speed and evolution of the waves, and they take a finite amount of energy 
out of the waves when they are accelerated and heated. This aspect of the 
calculation is particularly important in the case of He 

By starting the computation outside the corona, however, we do not 

deal with the problem of heavy ion abundance in the solar wind. We will 

assume that some process such as Coulomb friction pulls the heavier ions 

out of the corona ( Geiss, et al ., 1970; Joselyn and Holzer , 1978; Hollweg , 

1981h; Owockl , 1982) so that, at 10 R , the concentrations are such that sub- 

sequent motion produces observed values of n,/n at 1 AU* 

h p 

6) The electron density and bulk velocity are given by the requirements 
of quasi-neutrality and zero current flow. We neglect the electron inertial 
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terms and the thermal force (Braglnskil , 1965) so that the motion of the 
electron gas is governed only by pressure forces, the electric field, and 
Coulomb interactions vith the ions. The electron energy equation Includes 
a collisionless heat flux (Hollweg . 1976), q^ ■ 1.5 n^KT^v^^^jb. Ion heat 
fluxes are neglected. 

7) In the inertial frame, the center-of-mass velocity is taken to 
be radial, and the average magnetic field is carried out from the sun at 
this velocity. The tangent of the spiral angle of the average field is then 

W - ^ 

where B.^ +B^ "B^, fl_is the angular rotational velocity of the sun, 

9 r os 

and r is heliocentric distance. In this work, is defined in the corotating 
frame, but it follows from this radial assumption that the center-of-mass 
velocity in the inertial frame is simply 

8) Following the discussion of the previous section, the wave field 
is taken to consist of UKB, nondispersive Alfven waves propagating along 

-Y 

the average magnetic field. The wave power spectrum is a power law -k 
with a total variance of <6B^>. The propagation, dissipation and interaction 
with the solar wind plasma proceed according to the scenario and equations 
given previously, 
b) Equations 

We solve the steady-state, spherically-symmetric solar wind fluid 

equations in the corotating frame for the independent variables ^ohr* 

T , T , T. as functions of r. 
e p n 

The continuity, momentum, and energy equations for an ion species, i. 


N>|U> 
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under the above assumptions are 

n. V . r* ■ constant 
1 oir 
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(U) 


‘“’if? [<* ’olr*] ■ - J? 

+ + Pi ) J 

* ^“phr ” ”olr^ 

+ (X + W2) I Ky - »„1^) f (Py) (15, 


dT 


n.v . K , 

1 oir dr 


V . kT. . 
oir i dr 


dnj ic (Ti -Tj ) , 

- Q, + Z 3K.. ■ exp (-n..2) 

dr 1 ij m^^ + mj IJ 


* a * “*) hi ‘’ojr - 'oir)’ ' ("ij) 


(16) 


where is the solar mass, G Is the gravitational constant, and is the 
radial electric field. is taken equal to zero unless the waves are 
saturated. It is then determined by (from (3)) 

F d? <'’"Sr) ■ - 5 0l >r <\hr ' "ol?'' 


(17) 
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OF POOR QUALITY 

combined with the leading tern in equation (12) . 

Equation (15) is obtained by taking the component of the momentun 
equation parallel to the average magnetic field and multiplying by 
(r*b)“^ • /l + . The terms on the right-hand-side of (15) are, respec- 

tively, the centrifugal force in the corotating frame, the gravitational 
force, the thermal pressure force, the force from the radial electric 
potential, the non-dissipative Alfven wave pressure, the dissipative wave 
force, and the Coulomb friction. The right-hand-side of (16) contains 
the effects of heating due to wave dissipation, the energy exchanged 
with other species through Coulomb collisions, and the Joule heating from 
these collisions. The Coulomb terms are summed over all species j ^ 1, 
including the electrons. The symbols in the Coulomb terms are defined by 
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(18c) 


where the Coulomb logarithm has been taken to be InA ■ 21. Except for the 
terms involving Q^, equations (14) - (18) are equivalent to the system of 
equations written by McKenzie, et al . (1979); some typographical errors 


appearing in their paper have been corrected. The Coulomb collision terms 
are equivalent to Chose given by Burgers (1960). 
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The electron equations, subject to assumption (6), are 


■ i ’i“i 
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( 22 ) 


where here, the sunmiations are over all Ion species. The momentum equation 
(21) is used to determine the electric field. The first term of the energy 
equation (22) includes the divergence of the collisionless electron heat 
flux. 

With (21) inserted in (15) through the electric field, the equations 
(14) and (15) for each of the ion species, along with (22), form a set of 
five first order differential equations for the five variables: two radial 

velocities and three temperatures. Specification of the values of these 


quantities at the inner radius, taken to be 10 R^, will deteraine the solution. 
Since 10 Rg is beyond any critical points of these equations for the paraneters 
of interest, it is sufficient to sinply propagate the solution outward from 
this petition. The results of this calculation will be presented and dis- 


cussed in the next section. 
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IV Numerical Rtsulf and Dlicusslon 

In this section, we will consider the question of whether the resonant 
cyclotron interaction with high-frequency Alfvcn waves, in the "most favor- 
able case" scenario we have constructed, is capable of accelerating and heating 
the observed heavy ion populations in the solar wind to the observed conditions. 
In this investigation, we will deal primarily with alpha particles as the 
heavy ion species, since observations of this species are most extensive and 
e constraints are therefore strongest. 

We start with a "benchmark" solution for the helium case, using what we 
i^eel is a plausible set of initial parameter values at 10 R^. We take values 
fur the protons, electrons, and waves as suggested by the wave-driven high 
speed stream model of Hollweg (1978a) : 

V - 300 km/s T - 10^ ®K 

opr e 

T - 5 X 10** °K B • 1600 y 

P o 

n ■ 5000 cm“3 <6B^>/B^ • 0.2 (23a) 

P *0 

at 10 R . The waves will saturate when <6B^>/B^ - tj, and the spectral index 
a “0 

of Che wave power will be taken as f * 4. To test the ability of our most 
favorable case to preferentially accelerate and heat helium ions, we will 
start a typical concentration of helium at parameter values equal to Chose 
of the protons: 


n /n ■ 0.04 
o p 


V /v ■ T /T 
oar opr a p 


1 


at 10 Rg. 


(23b) 
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The velocities and temperatures from this benchmark case are shown in 
Figure 1. Figure la is a plot of the ion radial velocities » along with 
v^p^ + V^, where is the radial component of the Alfven velocity. 

The ion velocities, starting together, remain together until the waves 
saturate at 19 Rg. At this point, the resonant interaction begins to operate 
and the alpha particles are preferentially accelerated. However, it is clear 
from the figure that the calculated alpha velocity exceeds the proton velocity 
by only a fraction of the Alfven speed. Figure lb displays the electron, 
proton, and alpha temperatures and the behavior here is similar. The ion 
temperatures evolve together until the waves saturate. The resonant inter- 
action then heats both species strongly with preferential enhancement of 
T^ , but T^ / Tp at 1 AU is only 2.9, less than the "desired" value of 4. 

Since the benchmark parameters do not work, we investigate the effect 
of modifying these parameters. W*" *ind that it is quite difficult to change 
the character of the benchmark solution, in particular to produce 

Increasing the initial wave amplitude (or equivalently, decreasing the 
imposed saturation level) does not improve the situation and indicates 
further conflicts with observation. Since the waves accelerate the protons 
as well as the helium, the result of increased wave energy is to create 
faster solar wind speeds, with both i»;n species flowing faster and the dif- 
ferential speed essentially unchanged. This development is shown in Figure 2, 

where we plot the values of v and Av = v - v at 1 AU, alonz with 

opr ap oar opr ® 

4v^p at 60 Rg, as functions of the initial mean square wave amplitude, 

with all other parameters as in (23). The benchmark case is indi- 
cated by the vertical dashed line at * 5 x 10“® Gauss^. As <5|^> at 

10 Rg is increased from the benchmark value, the flow speeds are enhanced. 
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but the differential speeds are only slightly larger. For comparision» the 
Alfven speeds in the benchmark case are (60 Rg) « 103 km/s, (1 AU) 

• 29.8 km/s> and these values increase by a factor ^ V' 2 from the left edge 
of Figure 2 to the right. It is clear that .the inability of this model to 
produce larger differential speeds does not result from insufficient wave 
power. FurtheriL:)re, the strong correlation of with which is 

observed at 1 AU (Asbridge, et al .> 1976; Grunwaldt and Rosenbauer , 1978; 
Neugebauer, 1981a) is only vaguely approximated in these results. 

The effect of increased wave power on the ion temperatures is shown in 
Figure 3. Here, we plot the pruton and helium temperatures at 1 AU for the 
same cases used in Figure 2. To aid comparison of these results with obser*- 
vations, we have plotted the values against the proton velocity (see Figure 
2). The ion temperatures increase with v^ (and with the initial wave power) 
as expected, but decreases dramatically. This trend in is 

opposite to the observed positive correlation at 1 AU (Feynman , 1975; Bosqued , 
et al . , 1977; Neugebauer , 1981a; see Neugebauer , 1981b for compiled observa- 
tional results appropriate to Figures 2 and 3) . 

Investigation of other solutions shows that increasing the initiel values 
of T^, v^, or Y fails to yield differential speeds equal to the Alfven speed 
from this calculation. The effects of the wave dissipation are so pronounced 
that the initial velocities and temperatures of the ions are essentially 
irrelevant (if they are not so extreme as to make the dissipative interaction 
unnecessary), and the final relative ion parameters are determined by the 
dissipative interaction alone. 

The failure of the resonant cyclotron interaction to produce the 
observed differential speeds can be explained qualitatively by invoking this 
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OF poos 


dominance of the dissipative effects. We find that this interaction possesses 
a natural limit to the differential speeds, due to the simultaneous action 
of the waves on the protons* 

Assume that the flow d 3 mamics of all ion species are completely dominated, 
in some region, by the dissipative wave forces. Then, from (6)i the accelera- 
tion of a given ion species along the magnetic field is 


a^ 




I 


(24) 


Now consider the ratio of the wave acceleration of a heavy ion species to that 
of the protons, using the first term in (12) to give the relative heating rates. 
This gives 




-Y 


(25) 


Again, we note that the waves will not preferentially accelerate a heavy ion 

species with respect to the protons unless y>2. In addition, this mechanism 

will only preferentially accelerate an ion species up to a velocity such 

that - 1* Further wave acceleration will act equally on both species 

and will not increase the differential speed Av, . To estimate this critical 

hp 

differential speed, take m/k,. = v + V., and set the ratio (25) equal to 

11 op II A 

one. The critical differential speed is then 


Av, * = V. 
hp / 




(26) 


. i 
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and this is a limiting value for preferential acceleration of heavy ions by 

A 

cyclotron resonance with Alfven waves. For the benchmark case of alpha 
particles and y * bhls limiting differential speed is only « .29 

and it is clear from Figure la that this is too small to yield 
at any subsequent point inside 1 AU. The critical differential speed (26) 
is also a weak funtlon of Y> which explains why increasing Y in the cal- 
culation is not very effective in obtaining larger differential speeds. As 
Y «o for alpha particles, the critical differential speed approaches a 

maximum of V /2. 

A 

The critical differential speed in (26) was obtained by using the leading 
term in equation (12) which is equivalent to assuming very low parallel tem- 
peratures. This assumption is stretched somewhat, at the temperature peaks 
of the benchmark solution (since the calculation uses isotropic temperatures), 
so we return to the full heating rate in (12) and ask what effect the cor- 
rection term there has on the limiting speed (26). If y> 1» we find that, 
to lowest order, the hv^* will be larger (smaller) if c is positive (negative), 
where 

- « - 2"-p) 

This is an involved expression, but if the ion temperatures are proportional 
to the ion masses, the anisotropies are equal, and y^2, then the limiting 
differential speed will be lowered if the anisotropy is large enough: F>y/2- 

Since resonant ion-cyclotron heating generally increases F, we expect that 
warm plasma corrections to (26) will not result in significantly larger dif- 
ferential speeds. However, this point clearly needs further study. 
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Given the existence of this critical differential speed, it is easy 
to understand why the calculation is so insensitive to the initial helium 
parameters* Increasing v at 10 R„ simply means that the critical speed is 
reached sooner, resulting in less preferential heating and smaller helium 
temperatures. Any reasonable increases in at 10 are either eliminated 
by the Coulomb interaction before saturation or are simply swamped by the 
overwhelming wave heating. 

However, the critical differential speed in (26) has been derived by 
neglecting the action of thermal pressure forces compared to the dissipative 
acceleration. This is valid for most cases of interest, but it is possible 
to produce numerical results in which AVj^^ substantially exceeds the cri- 
tical differential speed through excessive pressure forces. This can be done 
by lowering the Initial heavy ion velocity, as seen in Figure 4. In this 
•’overshoot" case, the helium velocity at 10 Rg was reduced to one-half of 
the proton velocity, with the helium density adjusted to retain an abundance 
ratio of 4% near 1 AU. All other input parameters are given by (23). Ini- 
tially, Coulomb collisions act to equalize the ion velocities and produce 
some heating of the ions. However, at saturation the helium is still moving 
slower than the protons. The wave heating rate is a strongly increasing 
function of the difference between the wave phase speed and the bulk speed of 
the ion. Thus, after saturation the helium ions are heated and accelerated 
to an even greater extent than in the other calculations. The helium quickly 
overtakes the protons and the wave acceleration returns to normal, but by 
this point the helium has been so strongly heated that the pressure force 
causes the helium to overshoot the critical speed. In this case, however, 
the helium only reaches v * v 4* 2V./3 at 1 AU and still exhibits T < 4T 
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there. This result could be pushed further, but it seems unlikely that 
the helium would lag the protons by as much as 150 km/s at 10 R_, and the 
temperetures involved make the approximation for the heating rates of 
questionable accuracy. Still, the means by which heavy ions are propelled 
out of the corona are not well known and it is possible that further study 
of this overshoot phenomenon will prove valuable. 

Thus, with the possible exception of this overshoot case, we find that 
the resonant interaction with ion-cyclotron waves cannot accelerate or heat 
the solar wind helium population to the observed values with respect to the 
protons, even in the most favorable case. The fundamental difficulty is that 
the cyclotron resonance also acts on the solar wind protons, and the pre*- 
ferentlal interaction with helium is simply not preferential enough. 

Equation (26) also predicts that other heavy ions will find limiting 

differential speeds which depend on their mass-per-charge ratio. This is 

verified by the calculations displayed in Figures 5 and 6 for oxygen and iron 

ions, respectively. In these calculations, the proton and wave parameters 

were input as for the benchmark case (23), but the heavy ion abundance ratios 

were lowered to values comparable to those given by Bame, et al . (1979) and 

Mitchell, et al. (1981): n^/n_ » 5 x 10“^ and n„ /n_ = 5 x 10“^. The cal- 
c Op Fe p 

culation for each ion charge state was made separately and these abundances 
were taken for each case* At these densities, the oxygen and iron ions are 
essentially test particles, though in principle they are still allowed to 
affect the protons* We also caution the reader that the oxygen and iron 
results do not include helium, so the heavy ions are possibly given more 
energy than a full multi-ion calculation would yield. Thus, the velocities 
and temperatures shown in Figures 5 and 6 should perhaps be considered as 
upper limits. The results of these calculations at 60 and 1 AU are also 
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shown in Table 1, along with the benchmark results for He . 

The two most abundant oxygen ions. O’*"® and have mass-per-charge 

ratios only slightly larger than He"*^. Consequently, the velocities in 
Figure 5a are only slightly higher than the helium benchmark flow velocities. 
As expected from (26), the 0'**^ speeds are also marginally higher than those 
for 0'*'^, The oxygen temperature ratios are significantly lower at 1 AU 
than the T^/T > 16 which is often observed, and the ratios display the 

0 p 

same type of anticorrelation with proton speed as shown in Figure 3 for 

helium. These results are in conflict with the observations of oxygen in 

high speed solar wind (Ogllvie, et al ., 1980, 1982; Schmidt, et al .. 1980) 

which show v^ = v > v , T-/T “ 16, and the same positive correlation of 
0 a p 0 p 

T-/T with V as seen for helium. 

Op p 

The iron ions have a higher critical differential speed (26), and Indeed, 

they do reach speeds of v^ + in this calculation. However, the ratios 

T_ /T fall well below the mass proportionality level that is observed 
re p 

(Mitchell, et al .. 1981; Bochsler and Geis s. 1982), and the iron temperature 

ratios are also anticorrelated with v^. 

P 

We note, finally, that the temperature ratios with respect to helium, 
Tj^/T^, are more closely matched to the mass ratios than the equivalent values 
with respect to the protons (see Table 1). This is an intriguing result, 
indicating an intrinsic tendancy to mass-proportional heating by wave dissi- 
pation. It can be shown through an analysis similar to that used to derive 
the critical differential speed (26), that if the waves accelerate all the 
heavy ions to the same speed, the ions will be heated in proportion to their 
masses. This result is Independent of the expression for the dissipative 
heating, Q^, requiring only that the acceleration and heating be related 
by (24) and that the interaction take place quickly. In this picture, 
however, the protons, essentially carrying the waves, are not accelerated 
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to the same degree as the other Ions and so do not enter Into this mass- 
proportlonallty as we see from our model. Of course « in our model the ions 
are not all accelerated to the same speed, but rather to their critical 
speeds which depend on A^/Z^, so the mass-proportionality is not exact. 

However, ws can expect that a wave dissipation process which does accelerate 
all heavy ions to the same speed, as is observed (Ogilvie, et al », 1982) 
will necessarily produce the observed mass-proportional temperatures as well. 

In this study, we have concentrated on the interaction with left- 
polarized Alfven waves. A recent calculation by Marsch, et al . (1982c) has 
included effects of an interaction with right-polarized magnetosonic waves, 
and they find that their model can produce inside 30 Rg if 

they hypothesize enough power in the waves. 

There are a number of differences between their model and the work 
presented here. The model of Marsch. et al . (1982c) allows for anisotropic 
(bi-Maxwe Ilian) particle distributions, and uses dispersive waves with a 
phase speed in the plasma frame of 

a)/k||- V^|(l + m/Qp) [l + 4it (pj_ - P„)/Bj]| (27) 


where the minus (plus) sign refers to left- (right-) polarized waves and 
p^, p^^ are the proton pressures perpendicular to and parallel to the average 
field* They do not include a process equivalent to our saturation-and-cascade 
scenario^ so the required wave power at resonant frequencies is much larger 
than we have used. As a consequence of starting the calculation with large 
amounts of resonant wave power, the plasma is accelerated and heated almost 
immediately, over a distance of several hundredths of a solar radius rather 
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rather than the 1-10 interaction length we obtain. 

Since they use dispersive waves, Marsch, et al . (1982c) can only 

accelerate helium ions to differential speeds 4v <0.2 V, by the inter- 
op • A 

action with left-polarized ion-cyclotron waves, verifying statements to that 
effect in DH, McKenzie and Marsch (1982) and Section 2 of this paper. Sub- 
sequently, the plasma flows into regions of decreasing Alfven speed while 
the distribution functions evolve in the expanding heliospheric field. If 
the helium has been heated enough in the parallel direction, and if the 
Alfven speed decreases enough, the faster helium ions will move into cyclotron 
resonance with any right-polarized MHD waves that are present. These waves 
can further heat and accelerate the resonant ions. With equal intensities 
of right- and left-polarized waves, Marsch, et al . (1982c) produce results 
with 4v^p > well within 30 Rg, apparently demonstrating that cyclotron 
resonant interactions can explain the Hellos observations of helium. 

However, there may be a problem with this interpretation. Inserting 
the phase speeds (27) into (2), one obtains expressions for the parallel 
velocity needed by an ion in order to cyclotron resonate with waves of either 
polarization. These two expressions are never equal: the maximum velocity 

for left-polarized resonance is always less than the minimum velocity for 
right-polarized resonance. For a helium ion, this "resonance gap" is very 
large, being always greater than 2V^ + 4n(p^- P||)/Bg^^. At 30 Rg, 

for isotropic distributions, this gap can be greater than 600 km/s, correspond- 
ing to the thermal speed of a helium gas with a parallel temperature of 
T|l > 174 million degrees, (These numbers are even larger if Pj^> pj|, as would 
be expected from cyclotron heating by left-polarized waves.) A particle 
moving with a parallel velocity in this gap is not in resonance with either 
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mode, so it is clear that the resonant Interaction with the left-polarized 
waves cannot directly produce pairficiles in resonance with right-polarized 
waves. If particles traveling at the maximum velocity for left-polarized 
resonance subsequently move into a region of low enough Alfven speed, 
they can become resonant with the right-polarized waves, but Alfven speeds 
this low are typically found only beyond 1 AU. 

It is therefore puzzling that Marsch, et al . (1982c) find such a 
strong resonance with right-polarized waves Inside 30 R^, where the Alfven 
speed is so large. We suspect that their assumption of a Maxwellian distri- 
bution in parallel velocity is responsible. This assumption Insures that 
there are at least a small number of particles in resonance with the right- 
polarized waves at the start of the calculation. When these few particles 
are accelerated and heated by the right-polarized waves, they are interpreted 
as representing a widening of the entire velocity distribution and this 
widening is Interpreted as an increase in the number of resonant particles. 
Thus, the assumption that the distribution remains Maxwellian can artiflcally 
propagate particles across the resonance gap. Eventually, the entire dis- 
tribution may be propelled to high speeds by this numerically Induced process. 
This picture is only a conjecture, as yet without proof, but we feel this 
point needs further study. 
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V Summary and Conclusions 

We have investigated numerical solutions to the spherically symmetric, 
steady-state solar wind equations including acceleration and heating through 
a resonant interaction with left-polarized wa\'es. The equations for the 
resonant interaction model a scenario which we feel is a "most favorable 
case" for preferential acceleration and heating of heavy ions. Briefly, 
we assume dispersionless, parallel-propagating high-frequency Alfven waves 
which saturate when the magnetic fluctuations become comparable with the 
background field. The energy dissipated in this saturation process is 
assumed to cascade to resonant wavenumbers where it is picked up by the 
plasma through the quasi-llnear resonant cyclotron interaction. 

We find that, even in this most favorable case, the resonant cyclotron 

interaction is not capable of accelerating and heating the alpha particle 

population in the solar wind to the observed speeds and temperatures. We 

first point out that the quasi-linear interaction does not yield preferential 

heating or acceleration for heavy ions unless the wave power spectrum is 

steeper than k“^. Wave spectra this steep are rarely seen in the solar 

wind. Moreover, even if sufficiently steep spectra are assumed, we find 

that the resonant process in the case of helium, though acting preferentially 

on the helium, still accelerates protons to such an extent that the observed 

differential speeds, ®re not produced. We have succeeded in 

understanding this result analytically, in terms of a critical differential 

speed (Equation 26) for preferential acceleration. We also find that this 

most favorable case is unable to explain either the values of T /T , or Che 

o p 

T^/Tp - Vp correlation, which is frequently observed in high speed streams. 

The critical differential speed depends on the mass-per-charge ratio 
of the ion. Solar wind oxygen (primarily 0'*’^ and 0'*'^), with similar ratios 
to that of He"^, is affected in a similar manner to the helium. Iron ions 
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in the solar wind have higher mass'per-charge ratios, and therefore reach 

higher differential speeds than heliun or oxygen, if the wave spectre are 

steep enough to produce any preferential acceleration. This process is not 

able to produce the observed mass-proportional temperature ratios, T. /T ■ A 

n p 

at 1 AU. However, some tendency cowards mass-proportional heating is found 
among the heavy ions (excluding protons), yielding ^ 

In conclusion, we have shown that even a "most favorable case" model 
of the resonant cyclotron interaction with left-hand waves is not capable of 
producing the observed preferential acceleration and heating of heavy ions 


in the solar wind. 
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Fig. 1 


Fig. 2 


Fig. 3 


Fig. 4 


Fig. 5 


Fig. 6 


Figure Captions 

Benchmark case for helium, (a) Radial flow speeds of 
protons (chick line) and alpha particles (thin line) . 
Dashed line is the radial component of -I- V^. 

(b) Temperatures of protons (thick line ) , alpha parti- 
cles (thin line), and electrons (dashed line). 

Proton radial speed at 1 AU and differential speeds at 
60 Rg and 1 AU as functions of the initial wave power. 
Vertical dashed line indicates the location of the 
benchmark case. 

Ion temperatures and temperature ratio as functions of 
proton radial speed at 1 AU, obtained by varying the 
initial wave power. Vertical dashed line indicates 
the location of the benchmark case. 

Overshoot case for helium. Radial flow speeds (a) and 
temperatures (b) as in Fig. 1. 

Oxygen. Radial flow speeds (a) and temperatures (b) 
for protons, 0 , and 0 . 

Iron. Radial flow speeds (a) and temperatures (b) for 
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protons, Fe , Fe , and Fe 
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